Critical properties of the unconventional spin-Peierls system TiOBr 
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We have performed detailed x-ray scattering measurements on single crystals of the spin-Peierls 
compound TiOBr in order to study the critical properties of the transition between the incom- 
mensurate spin-Peierls state and the paramagnetic state at Tc2 ~ 48 K. We have determined a 
value of the critical exponent /3 which is consistent with the conventional 3D universality classes, in 
contrast with earlier results reported for TiOBr and TiOCl. Using a simple power law fit function 
we demonstrate that the asymptotic critical regime in TiOBr is quite narrow, and obtain a value 
of pasy = 0.32 ± 0.03 in the asymptotic limit. A power law fit function which includes the first 
order correction-to-scaling confluent singularity term can be used to account for data outside the 
asymptotic regime, yielding a more robust value of jSavg ~ 0.39 ± 0.05. We observe no evidence 
of commensurate fluctuations above Tci in TiOBr, unlike its isostructural sister compound TiOCl. 
In addition, we find that the incommensurate structure between Tci and Tc2 is shifted in Q-space 
relative to the commensurate structure below Tci. 

PACS numbers: 75.40.Cx, 78.70.Ck, 64.60.F- 



I. INTRODUCTION 

The physics of low dimensional quantum magnets is 
a very interesting and topical area of research. These 
systems exhibit a broad range of novel physical prop- 
erties and exotic quantum ground states, ranging from 
spin-Peierls and Haldane gap behavior in quasi-one- 
dimensional systems to high temperature superconduc- 
tivity in quasi- two-dimensional systems. Spin-Peierls 
systems in particular have undergone a recent renais- 
sance due to the discovery of two new, and somewhat 
unconventional, inorganic spin-Peierls compounds, the 
titanium oxyhalides TiOBr and TiOCl. 

Spin-Peierls behavior arises in quasi-one-dimensional 
(ID) systems which possess a combination of spin 1/2 
magnetic moments, short-range antiferromagnetic inter- 
actions, and strong magneto-elastic coupling. At low 
temperatures these systems undergo a spin-Peierls phase 
transition which is characterized by dimerization of the 
lattice and the development of a non-magnetic singlet 
ground state. As the presence of interchain interac- 
tions typically drives quasi-lD systems towards long- 
range magnetic order, only a select family of compounds 
with particularly strong spin-phonon coupling have been 
found to exhibit a spin-Peierls transition. The first ex- 
perimental realization of the spin-Peierls transition was 
discovered in organic charge transfer compounds, such 
as TTF-CuBDTi and MEM-(TCNQ)22. It was consid- 
erably later that spin-Peierls behavior was observed in 
CuGeOa, the first inorganic spin-Peierls compound (Tsp 
~ 14 K>g. 

TiOBr and TiOCl, the recently discovered inorganic 
spin-Peierls compounds based on Ti^+, are often re- 
ferred to as "unconventional" spin-Peierls systems. Un- 
like standard spin-Peierls systems, which undergo a 



single continuous phase transition at Tsp, the tita- 
nium oxyhalides undergo two successive phase transi- 
tions upon warming - a first-order transition at Tci be- 
tween commensurate and incommensurate spin-Peierls 
states, and a higher-order transition at Tc2 between 
the incommensurate spin-Peierls state and a disordered 
pseudogap state which extends up to T* ~ 135 K— . 
In addition, the titanium oxyhalides are unique among 
spin-Peierls compounds because of their unusually high 
transition temperatures (Tci/Tc2 ~ 65 K/92 K in 
TiOCli^i^ii^ and 27 K/48 K in TiOBr^iii) and sur- 
prisingly large singlet-triplet energy gap (Eg ~ 430 to 
440 K in TiOCliii^ and 150 K in TiOBpiS). The origins 
and properties of the incommensurate spin-Peierls state 
between Tci and Tc2 have attracted particular atten- 
tion to these compounds. 

Both TiOBr and TiOCl crystallize into a FeOCl-type 
structure, with an orthorhombic space group of Pmmn 
at room temperaturei^iii^. This structure consists of 
buckled Ti-0 bilayers in the a6-plane, which are sep- 
arated by double layers of either Br~ or Cl~ ions and 
stacked along the c-axis. The additional reduction of 
effective dimensionality within the Ti-0 bilayers, from 
two to one, is believed to arise from ordering of the Ti'^+ 
orbital degrees of freedom'^. Spin-Peierls dimerization 
occurs along the crystallographic 6-axis, and is driven 
by the direct exchange of Ti dxy orbitals^. High tem- 
perature susceptibility measurements suggest that both 
TiOBr and TiOCl can be fit to a S = 1/2 Heisenberg 
chain model, with nearest-neighbor exchange couplings 
of - 660 K and 360 K respectiveljs^ii^. 

In general, TiOBr has attracted considerably less at- 
tention than its sister compound TiOCl. This is largely 
due to the fact that it is more difhcult to grow single 
crystals of TiOBr than TiOCl, although it is also more 
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challenging to work with TiOBr because it is extremely 
hygroscopic (even more so than its isostructural coun- 
terpart) and has substantially lower transition temper- 
atures. All measurements to date seem to suggest that 
apart from a slight difference in energy scales, the two 
compounds behave almost identically - they have match- 
ing low temperature phase diagrams, and appear to be 
isostructural in each of the major phases. However, there 
are subtle differences between the lattice constants of the 
two materials, perhaps most notably in the ratio of the a 
and b lattice constants which describe the Ti-0 bilayers. 
While the room temperature lattice parameters reported 
for TiOBr arc a = 3.78 A, & = 3.49 A, and c = 8.53 AiS, 
those of TiOCl are a = 3.79 A, b = 3.38 A, and c = 8.03 
A". Thus, while the a-axis (i.e. interchain) spacing of 
the two compounds is almost identical, the 6-axis (i.e. 
intrachain) spacing is ~ 3 to 4% smaller in TiOCl. As 
a result, both the relative and absolute strength of the 
intrachain interactions are expected to be greater in Ti- 
OCl, giving rise to physical properties which are more 
strongly one-dimensional in nature. This appears to be 
reflected in both the magnitude of the nearest-neighbor 
exchange coupling and the size of the singlet-triplet en- 
ergy gap. 

The critical behavior of spin-Peierls systems is a sub- 
ject which has generated considerable interest and con- 
fusion over the past thirty years. Widely varying and 
often contradictory critical exponents have been re- 
ported for a number of spin-Peierls compounds. Early 
studies of CuGeOs, the first inorganic spin-Peierls sys- 
tem, led to conflicting claims of tricritica l^^'^^i^^ , three- 
dimensional^^'^"'^^, and mean-field-like^^ critical behav- 
ior. Subsequent work2^ has revealed that much of this 
discrepancy can be attributed to the narrowness of the 
asymptotic critical region, and it has been shown that 
the system belongs to a conventional 3D universality 
class. Similarly, while the organic compounds MEM- 
(TCNQ)2 and TTF-CuBDT were initially beheved to 
exhibit mean-field-like critical behavior— i^, more de- 
tailed measurements^ suggest that these compounds 
correspond to 3D universality classes as well. 

The critical behavior of the recently discovered Ti- 
based spin-Peierls compounds, TiOBr and TiOCl, has 
been less extensively studied but appears similarly enig- 
matic. To date there has been only one report of 
critical exponents for these systems, which placed the 
value of P for the continuous transition at Tc2 (between 
the disordered pseudogap state and the incommensu- 
rate spin-Peierls state) between 0.10 ± 0.04 and 0.26 
± 0.09 in TiOBr^. Similar values of /3 = 0.15^^ and 
P = 0.25^ can be extracted from other x-ray diffrac- 
tion measurements on TiOCl and TiOBr, although these 
results have not been explicitly quoted as critical ex- 
ponents by the authors. This range of /3 values sug- 
gests some form of either quasi-2D or tricritical behav- 
ior at Tc2, a result which is somewhat surprising for a 



phase transition which is known to be driven by inher- 
ently three-dimensional spin-phonon interactions. For 
two compounds which have already been shown to ex- 
hibit behavior which significantly differs from the stan- 
dard spin-Peierls scenario, this poses a very interesting 
question - is the critical behavior of these compounds 
truly distinct from that of CuGeOa and the organic spin- 
Peierls systems or is 3D universality one of the funda- 
mental properties of any spin-Peierls transition? The 
need for a detailed study of the critical behavior of 
these systems has been further underlined by a recent 
claim based on symmetry analysis^^ that the transition 
at T(72 is intrinsically first order in nature, despite the 
rather substantial weight of experimental evidence to the 
contrarj^i^iii^iSS^. 

In this paper we report the first detailed study of the 
critical properties of TiOBr. Using an approach which 
incorporates both standard power law fits and modified 
power law fits which include a first order correction-to- 
scaling term, we find that the transition at Tc2 can be 
well described by conventional 3D universality classes. 
We offer a qualitative comparison which suggests that 
TiOCl also corresponds to a similar universality class, 
and discuss the potential change in criticality which re- 
sults when TiOCl is doped with non-magnetic Sc ions. 
Our results indicate that the asymptotic critical regime 
in TiOBr is quite narrow, providing a natural explana- 
tion for the anomalously low values of (3 which have 
previously been reported in the literature. The ease 
with which our data can be fit to models of a contin- 
uous phase transition clearly contradicts arguments for 
the first order nature of the spin-Peierls transition at 
Tc2^. In addition, we offer evidence that commen- 
surate fluctuations, which inhabit the incommensurate 
spin-Peierls and pseudogap phases of TiOCl, are absent 
in TiOBr. We propose that this may reflect the weaker 
one-dimensionality of the system, and the greater im- 
portance of interchain interactions. Finally, we observe 
that the incommensurate and commensurate spin-Peierls 
states in TiOBr appear to be shifted in Q-space with re- 
spect to each other. This may be an indication that 
the structure of the incommensurate spin-Peierls state is 
more complex than previously assumed. 



II. EXPERIMENTAL DETAILS 

Single crystal samples of TiOBr were prepared using 
the chemical vapor transport method. A more detailed 
description of the crystal growth procedure has been re- 
ported elsewhere^.. The dimensions of the sample used 
in this experiment were approximately 4.0 x 3.0 x 0.5 
mm. The sample was mounted on the cold finger of a 
closed cycle refrigerator and aligned within a Huber four 
circle diffractometer. This sample environment provided 
a base temperature of ^ 7 K, with temperature stability 
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of better than ± 0.005 K. X-ray scattering measurements 
were performed using Cu-Ka radiation (A = 1.54 A) pro- 
duced by an 18 kW rotating anode x-ray source with a 
vertically focussing pyrolytic graphite monochromator. 



X-ray scattering scans were performed through sev- 
eral equivalent commensurate ordering wave vectors of 
the form (H, K-f 1/2, L), where superlattice Bragg peaks 
are expected to arise as the result of spin-Peierls dimer- 
ization. Due to the thickness of our sample and the rel- 
atively short penetration depth of the wavelength being 
employed 20 /im), these measurements were largely 
restricted to reflection geometry. Ultimately, our mea- 
surements focused around the (0, 0.5, 6) position in re- 
ciprocal space, as it was here that the intensity of the su- 
perlattice Bragg peaks was found to be strongest. Scans 
were carried out along the H, K, and L directions in re- 
ciprocal space in order to characterize all of the relevant 
commensurate and incommensurate superlattice peaks. 
These measurements were repeated at several tempera- 
tures within each of the four major regions of the TiOBr 
phase diagram - below Tci in the commensurate spin- 
Peierls phase, between Tci and Tc2 in the incommensu- 
rate spin-Peierls phase, between Tc2 and T* in the pseu- 
dogap phase, and above T* in the paramagnetic phase. 



Once the general features of the scattering at the su- 
perlattice peak positions had been well characterized, a 
more detailed study of the incommensurate scattering 
at (H ± 6h, K-I-1/2, L) was performed close to Tc2 in 
order to investigate the critical behavior of the system. 
In particular, a series of scans were taken through the 
incommensurate ordering wave vector at {-6h, 0.5, 6) to 
follow the temperature evolution of the incommensurate 
scattering intensity. To ensure that these measurements 
captured all of the relevant incommensurate scattering, 
two long scans were performed at every temperature - 
one scan in the H-direction to locate the center of the 
incommensurate satellite peak, followed by a second scan 
in the K-direction to obtain a measure of the integrated 
scattering intensity. Measurements were performed over 
a range of temperatures extending from ~ 11 K below 
the phase transition to ~ 3 K above the transition. The 
size of the temperature steps varied from as large as 0.25 
K (more than 6 K away from Tc2) to as small as 0.05 K 
(within 1 K of Tc2)- In order to check for the presence 
of thermal hysteresis at Tc2, two separate data sets were 
collected while warming and cooling through the transi- 
tion temperature. As no observable difference could be 
found between the measurements obtained on warming 
and cooling, both sets of measurements will be treated as 
a single data set in the analysis which follows in Section 
III.A. 
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FIG. 1: (Color online) Color contour maps of x-ray scatter- 
ing intensity as a function of temperature in TiOBr (left) and 
TiOCl (right). These maps are composed of H-scans of the 
form (H, 0.5, 6) and (H, 1.5, 1) respectively. A high temper- 
ature (150 K) background data set has been subtracted from 
the TiOCl data in order to eliminate higher order wavelength 
contamination. 



III. RESULTS AND DISCUSSION 

A. Critical Behavior of TiOBr 

Figure 1 shows two color contour plots, consisting of 
H-scans collected for TiOBr (left panel) and TiOCl (right 
panel). Each of these maps is composed of a series of H- 
scans performed through one of the superlattice peak 
positions - the (0, 0.5, 6) position in the case of TiOBr, 
and the (2, 1.5, 1) position in the case of TiOCl. Note 
that in the case of TiOCl a high temperature background 
data set, collected at T = 150 K, has been subtracted 
from each scan. This correction is necessary in order to 
eliminate weak A/2 contamination present in the inci- 
dent x-ray beam, which can give rise to Bragg scattering 
at the superlattice peak positions. The A/2 scattering at 
(2, 1.5, 1) in TiOCl is comparable in strength to the weak 
A scattering in the pseudogap phase (~ 0.1 counts/s), 
but it can be eliminated with a simple background sub- 
traction due to its lack of temperature dependence. No 
higher order wavelength contamination was observed at 
the (0, 0.5, 6) position in TiOBr, so a corresponding 
high temperature background subtraction has not been 
performed. 

The color contour maps in Fig. 1 offer a qualita- 
tive summary of the main features of the low tempera- 
ture phase diagram of the titanium oxyhalides. Below 
Tci there is a dimerized, commensurate two-fold su- 
perstructure, identified by a single superlattice peak at 
(H, K-I-1/2, L). Between Tci and Tc2 two incommen- 
surate satellite peaks appear at the (H ± 5h, K-I-1/2, L) 
positions, signifying the development of an incommen- 
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FIG. 2: (Color online) The temperature dependence of the 
integrated scattering intensity of the superlattice peaks ob- 
served at (0, 0.5, 6) [below Tci, in the commensurate SP 
phase] and {±5h, 0.5, 6) [above Tci, in the incommensurate 
SP phase] in TiOBr. The discontinuous drop in scattering 
intensity at Tci ~ 27 K is indicative of a first order phase 
transition, while the gradual decrease towards Tc2 is consis- 
tent with a continuous phase transition. 
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FIG. 3: (Color online) Fits to the integrated intensity of 
the incommensurate peak at (-&, 0.5, 6) as a function of 
temperature. These fits were performed using all data col- 
lected above Ti = 43 K, employing both a straight power 
law fit function (blue, jSasy ~ 0.32) and a power law includ- 
ing the first order correction-to-scaling confluent singularity 
term (red, /3avg = 0.39). In the inset the quality of these fits 
is shown on a logarithmic scale. 



surately modulated version of the dimerized structure. 
The primary component of the incommensurate modu- 
lation lies along the H-direction in reciprocal space, or 
the a- axis in real space, which means that it runs per- 
pendicular to the dimerized chains of Ti atoms in the 
spin-Peierls state. Note that additional commensurate 
scattering, which has been attributed to commensurate 
fluctuations^, is visible between Tci and T* in TiOCl 
but not in TiOBr. This result will be discussed in greater 
detail in section III.C. For the present time it is worth 
noting that, if anything, the physical picture for TiOBr 
is simpler than that of its isostructural counterpart. 

TiOBr is known to undergo two successive phase tran- 
sitions at low temperatures - a continuous transition into 
an incommensurate spin-Peierls state at Tc2 ~ 48 K, 
followed by a first order transition into a commensurate 
spin-Peierls state at Tci ~ 27 K. These two transitions 
are clearly illustrated in Fig. 2, which shows the temper- 
ature dependence of the integrated scattering intensity 
in the vicinity of the (0, 0.5, 6) superlattice peak po- 
sition. Below Tci the scattering intensity arises from 
the commensurate superlattice peak at (0, 0.5, 6), while 
between Tci and Tc2 it originates from the incommen- 
surate satellite peaks at (± Sh, 0.5, 6). All of the scat- 
tering intensities provided in Fig. 2 were obtained by 
integrating over H-scans of the form shown in Fig. 1. 
It should be noted that the intensity of the scattering 
which results from spin-Peierls dimerization is propor- 
tional to the square of the relevant atomic displacements, 
and hence can be directly related to the square of the or- 
der parameters associated with the phase transitions at 



Tci and Tc2. For this reason, the discontinuous jump 
in scattering intensity at Tci is evidence for a first order 
phase transition, while the gradual decrease towards Tc2 
is indicative of a continuous phase transition. Both of 
these results are consistent with the generally accepted 
nature of the transitions at Tci and Tc2- 

Since the intensity of the dimerization scattering is 
proportional to the square of the order parameter, it is 
possible to extract a value for the critical exponent /3 
from our x-ray diffraction data. As only the phase tran- 
sition at Tc2 is continuous, we will focus our attention 
on the critical behavior associated with the transition 
between the incommensurate spin-Peierls phase and the 
disordered pseudogap phase. In this case, the order pa- 
rameter for the phase transition will be related to the 
intensity of the incommensurate satellite peaks which 
appear at the (H ± Sh, K -I- 1/2, L) wave-vectors. The 
integrated scattering intensity of the {-Sh, 0.5, 6) incom- 
mensurate peak, which was found to be the strongest 
of the experimentally accessible incommensurate peaks, 
was carefully determined by performing a series of H and 
K-scans as described in Section II. The integrated scat- 
tering intensity was then fit to two different functional 
forms, a straight power law (Eq. 1) and a modified power 
law which includes the first order correction-to-scaling 
confluent singularity term (Eq. 2). 

I(t) = Aff^ + background (1) 
I{t) ^ At'^^ X (1 + Bf^) + background (2) 
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Here t = {T — Tc2)lTci is the reduced teniperature, 
while the value of the exponent A has been determined 
to be ^ Q_5 32,33 _ Representative fits to the integrated 
intensity data using Eqs. 1 and 2 are provided in Fig. 3. 
The choice of two fit functions was motivated by previ- 
ous studies on the critical behavior of CuGeOj^^ which 
showed that the asymptotic critical region was extremely 
narrow, extending less than 0.5 K (or ^ 4%) below T^p. 
In this region the length scale associated with fluctu- 
ations in the order parameter is greater than any of 
the other relevant length scales that occur in the sys- 
tem. Within the limit of the asymptotic critical regime 
it is expected that the order parameter should be ac- 
curately described by the simple power law provided in 
Eq. 1. Farther away from the transition temperature 
however, outside the asymptotic critical regime, other 
length scales in the system become significant and the 
behavior of the order parameter begins to deviate from 
this ideal form. It is this effect which was responsible for 
the wide variation between the initially reported critical 
exponents for CuGeOs, and we propose that the same 
holds true for TiOBr. In order to account for devia- 
tions from standard power law behavior farther away 
from Tc2, we can modify our original fit function to in- 
clude the first order correction-to-scaling confluent sin- 
gularity term, as shown in Eq. 2. It should be noted 
that this correction is only the first term in a larger se- 
ries expansion which describes the behavior of the order 
parameter beyond the limits of the asymptotic critical 
regime^. 

One of the factors that complicates the accurate deter- 
mination of critical exponents is the uncertainty inher- 
ent in simultaneously fitting both the critical exponent /3 
and the transition temperature Tc2 ■ As shown in Fig. 4, 
these two parameters are intimately related, and the op- 
timal value of /3 is strongly dependent on the appropriate 
choice of Tc2 ■ This situation can be even further compli- 
cated by the presence of critical fluctuations, which are 
typically found in the vicinity of a continuous phase tran- 
sition and can obscure the precise value of the transition 
temperature. Fortunately, as in the case of CuGeOg^, 
the critical fluctuations around in TiOBr are quite 
difhcult to observe. This is illustrated by the data pro- 
vided in Fig. 3, which clearly shows a sharp drop in 
scattering intensity at Tc2 with very little broadening 
or rounding near the transition. This makes TiOBr an 
excellent candidate for studies of critical behavior, and 
should allow for a particularly accurate determination of 
/3 and Tc2- 

The importance of the correction-to-scaling term is 
illustrated by the data provided in Fig. 5, which shows 
the effect of fitting range on the value of the critical 
exponents extracted from Eqs. 1 and 2. Here the fitting 
range is described in terms of the initial temperature, T^, 
which represents the start of the fitting range and the 
lower temperature bound for the data points included in 
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FIG. 4: (Color online) The dependence of the critical expo- 
nent /3 and goodness of fit parameter on the transition 
temperature Tc2. The bottom panel contains fit parame- 
ters obtained using a straight power law fit function (as de- 
scribed by Eq. 1), while the top panel shows parameters 
obtained from fits to a power law which includes the first or- 
der correction-to-scaling term (as described by Eq. 2). Both 
sets of fits were performed using data collected above Ti = 
43 K. 



the fit. Since the upper bound for each fitting range was 
fixed at a final temperature of T/ — 51 K, it follows that 
Ti can be taken as both a measure of the size of a given 
fitting range and an indication of how far a particular 
data set extends below the transition temperature Tc2- 
Thus higher values of T^ correspond to smaller fitting 
ranges, and describe fits to data which falls within (or 
at least close to) the asymptotic critical regime. 

The exponents determined from fits to the straight 
power law of Eq. 1 are shown by the black circles in 
Fig. 5. Note that these exponents follow a clear trend, 
gradually increasing as T^ rises and the fitting range is 
reduced. Closer to the transition temperature (i.e. above 
Ti 43 K) the values of (3 appear to level off, approach- 



ing an asymptotic limit of Pas 



0.32 ± 0.03. This 



behavior suggests that the asymptotic critical regime in 
TiOBr is quite narrow, extending less than 5 K (or ~ 9%) 
below T(72 ■ Within the critical region, where the critical 
behavior of the system can accurately be described by a 
standard power law, Eq. 1 provides consistent and re- 
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FIG. 5: (Color online) The value of the critical exponent 
/3 as a function of Ti, the start of the fitting range. The 
black circles represent values obtained from fits to a straight 
power law (Eq. 1), while the red squares represent values ob- 
tained from fits to a power law which includes the first order 
correction-to-scaling term (Eq. 2). The predicted values of 
/3 for the tricritical and conventional 3D universality classes 
have been provided as a guide to the eye. 



liable /3 values. Over larger fitting ranges, for data sets 
which extend beyond the critical region, deviations from 
simple power law behavior become increasingly signifi- 
cant and fits to Eq. 1 yield progressively worse values of 
f}. Hence Pasy represents the most accurate value of (3 
that can be obtained from Eq. 1, taken in the appropri- 
ate asymptotic limit close to Tc2- 

The critical exponents obtained from fits to the modi- 
fied power law of Eq. 2 are illustrated by the red squares 
in Fig. 5. In contrast to the results of the straight power 
law fits, the values of f3 obtained from fits to Eq. 2 follow 
a relatively constant trend, with slow variation around 
Pavg = 0.39 ± 0.05. Thus we find that the introduc- 
tion of the correction-to-scaling term compensates for 
deviations from standard power law behavior beyond the 
asymptotic critical regime, dramatically reducing the ef- 
fect of fitting range on the resulting critical exponents. 
Ideally, one might expect that sufficiently close to Tc2 
the values of /3 from the two fit functions should con- 
verge. In practice, the quality and reliability of fits per- 
formed when Ti is within 2 K of the phase transition 
are somewhat limited by the finite size of the data sets, 
which grow progressively smaller as T^ increases. How- 
ever, bearing this restriction in mind, it is very encour- 
aging to note that the values obtained for /3asy from Eq. 
1 and Pavg from Eq. 2 are fully consistent within the 
bounds of experimental uncertainty. 

From these results, we can conclude that the critical 
behavior of TiOBr can be described by one of the con- 
ventional 3D universality classes - 3D Ising (/3 = 0.326), 
3D XY (/3 0.345) or 3D Heisenberg (/? = 0.367). Both 



Pasy, the result obtained from straight power law fits 
performed in the appropriate limit, and Pavg, the result 
obtained from modified power law fits which included a 
correction-to-scaling term, are consistent with these val- 
ues. Due to the relatively large size of our error bars, it 
is difficult to narrow down the list of potential universal- 
ity classes to suggest a specific symmetry for the order 
parameter. One might only note that the best fits to the 
data are provided by the modified power law of Eq. 2, 
which leads to a critical exponent (3avg indicative of a 
3D universality class with continuous symmetry, either 
XY or Heisenberg. This would be consistent with pre- 
vious findings for conventional spin-Peierls systems such 
as CuGeOj^ and MEM-(TCNQ)2^. 

It should be noted that the success of our anal- 
ysis of the critical behavior in TiOBr provides ex- 
tremely compelling evidence that the phase transition 
at Tc2 is continuous, rather than first order, in na- 
ture. The fact that the integrated intensity of the 
superlattice peaks can be fit so well using a contin- 
uous model strongly contradicts the recent claims of 
Schonlcber et ai2^, while lending support to earlier find- 
ings based on magnetic susceptibihty'^'^, NMR"', spe- 
cific heali^'^, thermal expansion^, and x-ray diffraction 
measurement o^^'^^i^^'^" . There can be no confusing the 
results of the present study with the previously reported 
critical exponents of /3 = 0.10 to 0.26^^, 0.15^, and 
0.25ii for TiOBr and TiOCl. However, our examina- 
tion of the relationship between the fitting range em- 
ployed and the resulting value of /3 provides a simple and 
straightforward explanation for how these anomalously 
low critical exponents may have come about - they are 
the result of fits performed to measurements collected 
outside the asymptotic critical regime. 

The goodness of fit parameter, x^, is plotted as a 
function of fitting range in Fig. 6. Similar to Fig. 5, 
this plot serves to emphasize the importance of consid- 
ering correction-to-scaling terms when dealing with data 
sets which extend beyond the asymptotic critical regime. 
The most striking difference between the values of 
obtained from Eqs. 1 and 2 emerges over large fitting 
ranges (i.e. for low values of T^). For T^ less than 43 
K, as the fitting range grows larger one can observe a 
gradual monotonic increase in the values from the 
straight power law fits. Again, this reflects the fact that 
outside the asymptotic critical regime a simple power 
law form is no longer sufficient to describe the behavior 
of the order parameter. In contrast, the values from 
fits to the modified power law which incorporates the 
first order correction-to-scaling term remain small, and 
relatively constant, down to the lowest values of T^. For 
small fitting ranges (i.e. for T^ greater than 44 K), the 
quality of fit decreases markedly for both fit functions. 
This is simply a result of the reduced number of data 
points available for fits at higher Tj, which ultimately 
becomes the limiting factor for the accuracy of the fits 
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FIG. 6: (Color online) The goodness of fit parameter as 
a function of T^, the start of the fitting range. The black 
circles correspond to fits performed to a straight power law 
fit function (Eq. 1), while the red squares correspond to fits 
to a power law which includes the first order correction-to- 
scaling term (Eq. 2). 

in this regime. More generally, it is important to note 
that for any value of the goodness of fit is always bet- 
ter when the data is fit to Eq. 2 rather than Eq. 1. In 
other words, the correction-to-scaling fits consistently of- 
fer a better description of the data than a straight power 
law, regardless of the fitting range chosen. 

B. Implications for TiOCl 

Having completed our analysis of the critical proper- 
ties of TiOBr, it is natural that we next turn our atten- 
tion to what our results can suggest about the critical 
behavior of TiOCl. While we lack sufficient data to per- 
form the analysis necessary to precisely determine the 
value of f3 for TiOCl, we can clearly demonstrate that 
the critical exponents determined in the previous sec- 
tion provide a very convincing description of recently 
reported measurements on TiOCl^l. 

As illustrated by the colour contour maps provided in 
Fig. 1, the low temperature phase diagram of TiOCl 
is considerably more complicated than that of TiOBr. 
In particular, TiOCl is distinguished by the presence of 
commensurate fluctuations which persist from Tci up to 
T*, throughout both the incommensurate spin-Peierls 
and pseudogap phases. However, since the order pa- 
rameter associated with the transition at Tc2 is strictly 
related to the intensity of the incommensurate satellite 
peaks, rather than the sum of all dimerization scatter- 
ing, it is only the incommensurate scattering at (H ± 
Sh, K+1/2, L) which must be considered when studying 
the critical behavior of TiOCl. 

The upper panel of Fig. 7 shows the temperature 
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FIG. 7: (Color online) The temperature dependence of the 
integrated scattering intensity at the (2 ± Sh, 1-5, 1) su- 
perlattice peak position in TiOCl (top) and Tio.ggSco.oiOCl 
(bottom). These intensities have been obtained by integrat- 
ing over H-scans of the form shown in the colour contour 
maps of figure 1. The TiOCl data set can be well described 
by a straight power law fit which is consistent with conven- 
tional 3D universality (/3 = 0.34, shown in blue), while the 
Tio.ggSco.oiOCl data set is much better fit by a power law 
function consistent with mean-field-like behavior (/3 = 0.50, 
shown in red). 



dependence of the integrated scattering intensity of the 
(2 ± 6h, 1.5, 1) incommensurate satellite peaks in Ti- 
OCl. These intensities were obtained by integrating 
over small Q-ranges around the incommensurate order- 
ing wave vectors in H-scans similar to those used to con- 
struct the colour contour maps of Fig. 1. The data 
were fit to a standard power law, as described by Eq. 
1, although due to the low density of points near Tc2 
the value of the transition temperature was fixed at 92 
K. This fixed value of Tc2 is consistent with previously 
reported susceptibility^ii, heat capacity^i^, NMR^, and 
x-ray diffraction measurements^'^^'^ for TiOCl. The re- 
sulting line of best fit (indicated by the solid blue line in 
the upper panel of Fig. 7) corresponds to a critical ex- 
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ponent of P — 0.34 ± 0.04. It should be noted that this 
value of (3 is consistent with the theoretical predictions 
of the 3D Ising, 3D XY and 3D Heisenberg models, indi- 
cating that the transition in TiOCl can be well described 
by any of the conventional 3D universality classes. Fur- 
thermore, there is excellent agreement between the value 
of /3 = 0.34 ± 0.04 for TiOCl and the value of f3asy = 
0.32 ± 0.03 obtained earlier from power law fits to the 
TiOBr data. In fact, this relatively simple estimate for 
/3 also agrees quite well with the value of (iavg — 0.39 ± 
0.05 obtained for TiOBr using the modified power law 
of Eq. 2, in spite of the fact that the large fitting range 
used for the TiOCl data set almost certainly extends well 
beyond the asymptotical critical regime. 

Unfortunately, due to the small size and low point 
density of the TiOCl data set, it is difficult to extract 
any meaningful parameters from fits performed using the 
modified power law described by Eq. 2. On a more qual- 
itative level however, it is worth noting that the function 
in Eq. 2 can fit the TiOCl data very well using a fixed 
transition temperature of Tc2 = 92 K and the critical 
exponent of Pavg — 0.39 determined for TiOBr. By eye, 
such a fit is almost indistinguishable from the solid blue 
curve which represents the best fit line for Eq. 1, hence 
it has been omitted from the upper panel of Fig. 7. 

A very interesting comparison can also be made using 
previous x-ray data collected on Tip.ggSco.oiOCl, a com- 
pound in which Sc-doping has been shown to suppress 
commensurate spin-Peierls order to below at least T ~ 
7 K— . The temperature dependence of the integrated 
intensity of the incommensurate scattering at (2 ± Sh, 
1.5, 1) in Tio.ggSco.oiOCl is provided in the bottom panel 
of Fig. 7. If we attempt to fit this data using the simple 
power law fit function described by Eq. 1 we find that 
the critical behavior of Sc-doped TiOCl is clearly differ- 
ent from that of either pure TiOCl or TiOBr. The best 
power law fit to the Tig.ggSco.oiOCl data (illustrated by 
the solid red line in the bottom panel of Fig. 7) corre- 
sponds to a critical exponent oi (3 = 0.50 ± 0.07 and a 
transition temperature of Tc2 = 109 ± 6 K. This value of 
/3 is consistent with the exponent predicted by mean field 
theory (/3 = 0.5), and is incompatible with the theoreti- 
cal predictions for any of the conventional 3D universal- 
ity classes. In addition, the value of /3 for Tig.ggSco.oiOCl 
is dramatically different from the critical exponents that 
we determined for TiOCl and TiOBr. This difference in 
critical behavior is evident even from the raw data pro- 
vided in Fig. 7. The intensity of the incommensurate 
scattering in Tio.ggSco.oiOCl appears to grow linearly 
with decreasing temperature, while the scattering in Ti- 
OCl clearly follows a sublinear power law trend. 

It is also instructive to compare these results to pre- 
vious studies on the critical properties of CuGeOg^i^. 
It has been well established that in both pure CuGeOa 
and CuGeOa doped with ions of similar ionic radius the 
critical behavior of the spin-Peierls transition can be de- 



scribed by one of the conventional 3D universality classes 
with continuous symmetry (such as 3D XY or 3D Heisen- 
berg). However, when CuGeOs is doped with ions of 
larger ionic radius, by substituting Cd^"*" for Cu^"*" or 
Ga''+ for Ge'*+ for example, the observed universality 
class is found to change towards mean-field-like behav- 
ior. This effect appears to have little or no concentration 
dependence, and is believed to be related to the presence 
of local strain fields induced by the larger dopant ions. 
The significance of local strain effects on the behavior of 
the order parameter is reasonably intuitive given both 
the magneto-elastic nature of the spin-Peierls transition 
and the underlying importance of lattice properties in 
the spin-Peierls materials. Hence, if the dopant ion is 
smaller or of comparable size to the ion it is replacing (i.e. 
if Aion is less than roughly -1-3% of the ionic radius) the 
universality class of the transition remains unchanged. 
Alternatively, for large, positive changes in ionic radius 
(i.e. Aion is on the order of -1-9% or more) there ap- 
pears to be a distinct change in critical properties. The 
difference in ionic radius between Sc^+ and Ti'^+ cor- 
responds to Aion = +11%. Thus, the change in critical 
behavior we observe between TiOCl and Tig.ggSco.oiOCl 
appears fully consistent with lessons learned from stud- 
ies of doped CuGeOg^. Although further study is still 
required to provide a more detailed understanding of this 
phenomenon, on the basis of these results it is certainly 
tempting to conclude that doping-induced changes in 
criticality are a fundamental property associated with 
the spin-Peierls transition. 



A second significant feature which arises from the com- 
parison of the TiOCl and Tig.ggSco.oiOCl data sets is the 
change in transition temperature, Tc2, which is induced 
by Sc-doping. As shown by Fig. 7, the value of Tc2 
is considerably higher in Tio.ggSco.oiOCl than in TiOCl, 
with the fitted value of Tc2 = 109 ± 6 K for the doped 
sample lying well above the accepted value of 92 ± 1 
K for the pure sample. This result is quite surprising, 
particularly when one considers that doping has been 
universally found to suppress the spin-Peierls transition 
temperature in CuGeOs. One possible explanation for 
this effect is that since no commensurate fluctuations are 
observed below T* ~ 135 K in Tio.ggSco.oiOCl, the T ~ 
109 K temperature scale actually corresponds to the ap- 
propriately lowered equivalent of T* from TiOCl. This 
argument is supported by the fact that the incommen- 
surate scattering in Tio.ggSco.oiOCl is very broad, im- 
plying that the length scale for correlations in the doped 
compound is even shorter than that of the commensurate 
fluctuations in pure TiOCl. In addition, the ^ 20 % drop 
between the 135 K and 109 K temperature scales is quite 
consistent with the 15 to 45 % suppression of Tsp per 1 
% concentration of dopants which has been reported in 
the literature for doped CuGe03^i2ii38^^^^^_ 
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C. Details of the Incommensurate Structure of 
TiOBr 

As discussed earlier, there are several respects in which 
the low temperature scattering of TiOBr differs from 
that of its isostructural counterpart TiOCl. The chief 
feature which distinguishes these two compounds from 
each other is the presence of commensurate fluctuations, 
which are clearly evident in TiOCP^ but appear to be 
absent in TiOBr. In TiOCl, these commensurate fluc- 
tuations can be observed throughout both the incom- 
mensurate spin-Peierls phase (from Tci to Tc2) and 
the pseudogap phase (from Tc2 to T*), as shown in the 
right panel of Fig. 1. The commensurate scattering is 
not resolution-limited above T^i, indicating that it does 
not correspond to static long-range order. However, the 
width of the scattering implies the presence of relatively 
long correlation lengths, on the order of ^ 100 A or more. 
Interestingly, between Tci and Tc2 the commensurate 
fluctuations appear to both coexist with and compete 
against the incommensurate scattering at (H ± iJ//, K -|- 
1/2, L). The commensurate fluctuations also appear to 
be shifted in Q-space with respect to both the accom- 
panying incommensurate scattering and the long-range 
ordered commensurate state below Tci, with observed 
peak positions displaced along both the H and K direc- 
tions. 

In TiOBr, we find no indication of commensurate fluc- 
tuations in either the incommensurate spin-Peierls state 
(between Tci and Tc2) or the pseudogap state (be- 
tween Tc2 and T*). This was the case at each of the 
7 different commensurate ordering wave-vectors investi- 
gated in this study. One possible explanation for the ab- 
sence of commensurate fluctuations in TiOBr may arise 
from the subtle differences between the structures of Ti- 
OCl and TiOBr. While the distance between neigh- 
boring dimerized chains of Ti atoms (i.e. the inter- 
chain spacing, related to the a lattice parameter) is 
almost identical in both two compounds, the distance 
between neighbouring Ti atoms on the same dimerized 
chain (i.e. the intrachain spacing, given by the b lat- 
tice parameter) is approximately 4 % larger in TiOBr. 
As a result of this change in lattice parameters, TiOBr 
appears to be less one-dimensional than its isostruc- 
tural counterpart. The difference is clearly evident from 
the magnetic susceptibility measurements reported for 
the two compounds^ i^°'^^ , both in terms of the magni- 
tude of the nearest-neighbor exchange coupling (which 
is estimated to be 45 % smaller in TiOBr) and the 
observed deviations from the idealized Bonner-Fisher 
curve (which suggest signiflcantly stronger interchain 
couphngs in TiOBr). Doping studies on TiOCl5i indi- 
cate that the relative strengths of the intrachain and in- 
terchain interactions are intimately connected with the 
stability of the incommensurate structure. Thus, if the 
interchain coupling is stronger with respect to the intra- 



chain coupling in TiOBr than it is in TiOCl then it could 
easily lead to suppression of commensurate fluctuations 
(as in Tio.ggSco.oiOCl), stabilization of the incommensu- 
rate structure, and a loss of phase coexistence between 
Tci and Tc2- 

It should be emphasized that the lack of observable 
commensurate fluctuations in the present study does not 
necessarily imply that such scattering is completely ab- 
sent in TiOBr. In fact, the difficulty in detecting such 
a weak, broad feature - particularly one which may be 
slightly shifted in Q-space - is clearly underscored by 
the number of x-ray diffraction studies which previously 
failed to observe such commensurate fluctuations in Ti- 
OCl. It should also be noted that the (2, 1.5, 1) wave- 
vector, where the commensurate scattering was found 
to be strongest in TiOCl, was experimentally inaccessi- 
ble in this study. Thus, there remains a possibility that 
stronger commensurate fluctuations could still be found 
in TiOBr at other, as of yet unstudied, positions in re- 
ciprocal space. Further work is needed to provide a more 
thorough understanding of this intriguing difference be- 
tween the TiOX spin-Peierls systems. 

Our measurements in the incommensurate spin-Peierls 
phase also reveal an interesting property of the in- 
commensurate modulation in TiOBr. In addition to 
the primary component of the incommensurate mod- 
ulation wave-vector which runs along the H direction 
(perpendicular to the direction of the dimerized chains), 
a smaller secondary component of the modulation has 
been reported along the K direction (parallel to the 
dimerized chains) in both TiOCl and TiOBr. The mag- 
nitude of the K-component of the incommensurate mod- 
ulation, 5k, increases monotonically between Tci and 
Tc2, reaching a maximum value of ~ 0.015 r.l.u. in 
both compoundsSSi^. Even though the experimental 
resolution of our present study is too broad to fully re- 
solve the distinct incommensurate satellite peaks which 
should arise from the K-component of the incommensu- 
rate modulation, we would anticipate that the effects 
of this secondary incommensuration should be mani- 
fested as a significant broadening of the superlattice 
peaks along the K-direction. However, as the represen- 
tative K-scans provided in the upper panel of Fig. 8 
clearly demonstrate, there is no observable change in the 
width of the {-Sh, 0.5, 6) superlattice peak between the 
commensurate and incommensurate spin-Peierls phases. 
What these K-scans do reveal is an obvious shift in the 
position of the incommensurate satellite peaks above 
Tci relative to the commensurate superlattice peaks be- 
low Tci- The temperature dependence of this shift in 
Q-space, which we will call Sk, is illustrated in the bot- 
tom panel of Fig. 8. 

Note that superlattice reflections are observed at the 
(± Sh, 0.5 + 5k, 6) positions, but not at (± 5h, 0.5 - 5k, 
6). In this respect, our K-shift clearly differs from the 
previously reported K-incommensuration, which when 
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FIG. 8: (Color online) Representative K-scans performed 
through the {-5h, 0.5, 6) superlattice peak in TiOBr (top). 
The black circles correspond to data collected within the com- 
mensurate spin-Peierls state (T — 20 K, 5h ~ 0), while the 
red squares correspond to data collected within the incom- 
mensurate spin-Peierls state (T = 45 K, Sh ~ 0.08). For 
illustrative purposes the intensity of the T = 45 K data set 
has been scaled by a factor of four. The solid lines repre- 
sent fits to a simple two-Gaussian function. Note that the 
peaks in the incommensurate phase are shifted in Q-space 
relative to those of the commensurate phase. The temper- 
ature dependence of the shift, Sk, is shown in the bottom 
panel, along with the H-component of the incommensurate 
modulation wave vector, Sh- 



coupled with the incommensuration along H should pro- 
duce a set of 4 incommensurate peaks at (± 6h, 0.5 
± Sk, 6). Furthermore, there are indications that this 
K-shift exhibits some form of systematic Q-dependence. 
At the (0, 0.5, 6), (0, -0.5, 6), and (1, 0.5, 6) positions 
in reciprocal space we find the incommensurate peaks 
shift to higher Q positions (i.e. to K -I- Sk), while 
for (0, 0.5, 5) and (1, 0.5, 5) we find the peaks shift 
to lower Q (i.e. to K - Sk)- At the (0, 1.5, 6) wave- 



vector no discernable K-shift could be observed. Based 
on alignment scans performed through several nearby 
structural Bragg peaks we can conclude that this shift 
in Q-space is not simply the result of a change in the 
lattice constants of the material. Interestingly, the shift 
in Q-space appears to have a different temperature de- 
pendence than the H-component of the incommensurate 
modulation wave-vector, as shown in the bottom panel 
of Fig. 8. While Sh and Sk both increase steadily with 
temperature above Tci , Sh approaches saturation by ^ 
35 K, while Sk continues to grow until the system reaches 
Tc2- If the observed K-shift was simply equivalent to the 
K-component of the incommensurate modulation wave- 
vector it seems surprising that it would not follow the 
same temperature evolution as SH. 

This effect has also been observed in recent x-ray 
measurements performed by Sasaki et ali^, which show 
that the incommensurate scattering at (± Sh, 2.5, 0) 
is shifted in the K-direction (towards lower Q) with re- 
spect to the commensurate scattering at (0, 2.5, 0). Fur- 
thermore, the magnitude and temperature dependence of 
this shift in Q-space are almost identical to that reported 
here. Curiously, in spite of the absence of additional 
satellite peaks at (± Sh, 2.5 -I- Sk, 0), Sasaki et al. treat 
this behavior as equivalent to the K-incommensuration 
reported by van Smaalen et alfii. This assumption is 
not completely unjustified, particularly since recent x- 
ray diffraction measurements on TiOCl^^ have shown 
that in some cases there can be dramatic differences be- 
tween the observed scattering intensities of neighbouring 
pairs of incommensurate satellite peaks. These intensity 
variations have been attributed to changes in structure 
factor resulting from phase shifts between the modula- 
tions of neighboring chains of Ti atoms in the unit cell. 
In TiOCl these intensity differences have been found to 
be up to one to two orders of magnitude in sizc^^. If 
the scattering at the (± Sh, 0.5 - Sk, 6) position in 
TiOBr is suppressed by a similar factor then it could be 
too weak to be detectable with our current experimental 
configuration. However, while it is possible that weak 
superlattice scattering could be obscured by experimen- 
tal background in the present study, it is less likely that 
this same problem would also have affected the results 
of Sasaki et al. due to the much higher signal-to-noise 
ratio available at synchrotron x-ray sources. 

An alternative possibility, which naturally arises from 
the parallels between TiOCl and TiOBr, is that the K- 
shift between the commensurate and incommensurate 
scattering in TiOBr is an analog of the Q-shift be- 
tween the commensurate fluctuations and long-range or- 
dered commensurate/incommensurate scattering in Ti- 
OCl. The magnitude of the Q-shift reported for the 
commensurate fluctuations in TiOCP^ is very similar to 
the size of the K-shift we observe in TiOBr, and it is clear 
that both of these features are intimately connected with 
the flrst-order phase transition at Tci. Thus, the shift 
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in Q-space between the commensurate and incommen- 
surate scattering in TiOBr could represent a distinctly 
different effect from the K-component of the incommen- 
surate modulation. In this case, the K-shift could reflect 
changes in the balance of interchain and intrachain inter- 
actions at Tci. Hence, below and above Tci the ratio of 
relevant interaction strengths favors commensurate and 
incommensurate modulations of the dimerized structure 
with two slightly different wave-vectors. It is difficult to 
comment on the specific implications of these Q-shifts 
without further, and more thorough, study of additional 
superlattice peak positions. However, this result may 
provide an indication that the structure of the incom- 
mensurate spin-Peierls state of TiOBr may be more com- 
plex than previously assumed. 

IV. CONCLUSIONS 

In conclusion, we have performed the first detailed 
study of the critical properties of TiOBr. While the 
spin-Peierls behavior of this system appears to be un- 
conventional in a number of respects, the critical prop- 
erties of the system can be very well described by conven- 
tional 3D universality classes. This finding is consistent 
with previous results obtained from studies performed on 
other spin-Peierls systems such as MEM-{TCNQ)-2^ and 
CuGeOa^^. By analyzing our data using a simple power 
law fit function we determined that the critical exponent 
associated with the transition at Tc2 approaches a value 
of Pasy = 0.32 ± 0.03 in the appropriate asymptotic crit- 
ical regime. The strong dependence of (3 on the size of 
the fitting range suggests that the asymptotic critical 
region in TiOBr is quite narrow, as in CuGeOg^. This 
provides an explanation for the apparently contradictory 
critical exponents that have previously been reported for 
this compound. The narrowness of the critical regime 
can be accounted for by employing a modified power law 
fit function which includes the first order correction-to- 
scaling confluent singularity term. Using this approach 
we obtain a value of (3avg = 0.39 ± 0.05, a result which 
is consistent with both the 3D XY and 3D Heisenberg 
universality classes. 



Based on a simple, primarily qualitative comparison, 
we argue that TiOCl, the isostructural cousin of TiOBr, 
can also be well described by one of the conventional 3D 
universality classes. More detailed studies would still be 
required in order to offer a more quantitative descrip- 
tion of the critical behavior of TiOCl, or to further nar- 
row down the list of potential universality classes for the 
transition. However, this result appears to add growing 
weight to the notion that the spin-Peierls transition is 
by nature a conventional three-dimensional phase tran- 
sition, whether it is realized in organic or inorganic sys- 
tems or even the "unconventional" titanium oxyhalides. 
Applying the same approach to measurements performed 
on Tio.ggSco.oiOCl we flnd evidence that Sc-doping in- 
duces a fundamental change in the critical properties 
of the transition at Tc2- We suggest that this effect 
is likely caused by the presence of local lattice strains 
which arise due to the signiflcantly larger ionic radius of 
the Sc^^ dopant ions. 

In addition, we describe several interesting features 
of the incommensurate structure of TiOBr between Tci 
and Tc2 ■ Our measurements reveal no evidence for com- 
mensurate fluctuations of the type observed in the in- 
commensurate spin-Peierls and pseudogap phases of Ti- 
OCl. This result may reflect the weaker one-dimensional 
properties of TiOBr, and in particular the substantially 
smaller size of the intrachain interaction. We also ob- 
serve that the commensurate and incommensurate struc- 
tures are shifted in Q-space with respect to one another, 
with a displacement along the K-direction that steadily 
grows between Tci and Tc2- We hope that these mea- 
surements will help to guide and inform future theoret- 
ical and experimental studies of these novel spin-Peierls 
systems. 
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